Effects of disorder on the non-zero temperature Mott Transition 
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The physics of the metal-insulator coexistence region near the non-zero temperature Mott tran- 
sition is investigated in presence of weak disorder. We demonstrate that disorder reduces the tem- 
perature extent and the general size of the coexistence region, consistent with recent experiments 
on several Mott systems. We also discuss the qualitative scenario for the disorder-modified Mott 
transition, and present simple scaling arguments that reveal the similarities to, and the differences 
from, the clean limit. 
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I. INTRODUCTION AND MOTIVATION 

The physics of the metal-insulator transition has con- 
tinued to attract considerable interest in recent years. 
Substantial progress has been achieved in understand- 
ing the behavior near the interaction-driven transition, 
where Dynamical Mean Field Theory 1 (DMFT) has been 
very successful in explaining the behavior of several 
classes of materials ranging from transition metal oxides 
such as V2O3 to organic Mott systems. This approach 
has been especially useful in describing the non-zero tem- 
perature behavior in the paramagnetic coexistence region 
between the metal and the insulator. In this regime, the 
two phases compete, and the resulting behavior emerges 
as a compromise between the energy gain to form coher- 
ent quasiparticles, and the larger entropy inherent to the 
incoherent insulating solution. There is not actual two- 
phase coexistence (as in conventional first order thermo- 
dynamic phase transitions) in this region. Rather, it is 
a region of parameters in which two local minima of the 
free energy coexist. 

So far, most theoretical work has concentrated on clean 
systems, although several experimental studies indicate 
that effects of disorder are particularly important pre- 
cisely in this coexistence regime. Measurements per- 
formed in compounds such as NiSSe mixtures 2 * 3 *^ and k- 
organics 5 '- indicate that the presence of disorder pushes 
down the critical temperature end point of the metal and 
insulator coexistence region. In particular, experiments 
performed on a N1S2 compound, which has much weaker 
disorder, show that the Mott transition occurs at 150K, 2 
with an external applied pressure of 3 GPa, while in the 
substituted NiS2-xSe ;E compound it is seen only below 
100K. 4 It is important to notice that applying an external 
pressure to these compounds is equivalent to substitut- 
ing S by Se, which might suggest that the results above 
would be in conflict. A speculation was made that the 
reduction in the transition temperature would be due to 
the local randomness introduced with Se substitution. 3 

We address the theoretical issues from the perspective 
of the Hubbard model. It is not a priori obvious what 
should be the effect of disorder on the size and the tem- 



perature range of the coexistence region. On the one 
hand, disorder tends to broaden the Hubbard bands and 
thus larger interaction is needed to open a Mott Hubbard 
gap. This may lead to a larger overall energy scale, which 
could stabilize the coexistence region. On the other hand, 
disorder generally leads to spatial fluctuations in all lo- 
cal quantities, an effect that could smear or decrease the 
jump at any first order phase transition, and thus reduce 
the coexistence energy scale. These considerations indi- 
cate that careful theoretical work is called for, which can 
address the interplay of interactions and disorder near 
the Mott metal-insulator transition. 

A formalism that describes the effects of disorder 
within a DMFT approach was outlined some time ago^i 
but a very limited number of calculations were explicitly 
carried out within this framework. More recently, the 
approach was reexamined to investigate strong correla- 
tion effects on disorder screening^ and the related tem- 
perature dependence of transport in the metallic phase, 9 
These results shed light on several puzzling phenomena 
observed in experiments on two dimensional electron sys- 
tems, but did not provide a description of the physics rel- 
evant to the coexistence region at non-zero temperature. 

In this paper we examine the phase diagram for the 
Mott transition in the presence of moderate disorder at 
non-zero temperature within the DMFT approach^ We 
present results describing the evolution of the coexistence 
region, showing that disorder generally reduces its size, in 
agreement with experiments. Our results give a physical 
picture that describes the gradual destruction of quasi- 
particles as the Mott insulator is approached, and estab- 
lish the qualitative modification of the critical behavior 
resulting from the presence of disorder. 

Our findings are valid in the regime of strong corre- 
lations but weak to moderate disorder, where Anderson 
localization effects, which are neglected in our theory, 
can be safely ignored. The latter have been included in 
earlier zero temperature DMFT-based strong correlation 
calculations In particular, we mention that our low- 
est temperature results are consistent with the T = 
result at weak disorder of Byczuk et aljii but give the 
temperature dependence of the metal-insulator coexis- 
tence region. 
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II. NON-ZERO TEMPERATURE DMFT FOR 
DISORDERED ELECTRONS 

We consider a half-filled Hubbard model in the pres- 
ence of random site energies, as given by the hamiltonian 



H 



<ij><7 



e l n i 



(1) 



Here cj (cj CT ) creates (destroys) a conduction electron 
with spin a on site i, iiity — Cj_Ci CT is the particle number 
operator, i is the hopping amplitude, and U is the on- 
site repulsion. The random site energies £j are assumed 
to have a uniform distribution of width W . 

Within DMFT for disordered electrons^ a quasiparti- 
cle is characterized by a local but site-dependent^ self- 
energy function £j(w) = £(cl>,£;). To calculate these 
self-energies, the problem is mapped onto an ensemble 
of Anderson impurity problems-' embedded in a self- 
consistently calculated conduction bath. In this ap- 
proach, only quantitative details of the solution depend 
on the details of the electronic band structure; in the fol- 
lowing we concentrate on a semi-circular model density of 
states. In this particular case, the hybridization function 
is given by 



A(w) = t 2 G{uj) 



(2) 



and the average local Green's function, G(u>), is obtained 
by imposing the following self-consistent condition 



1 



A(w) - S;M 



(3) 



where (...) indicates the arithmetic average over the dis- 
tribution of £,-. 

To solve the single-impurity problems at non-zero tem- 
perature for different site energies, we mostly used the 
iterated perturbation theory (IPT) method of Kajueter 
and Kotliarpi&ii However, to check the accuracy of the 
results, in several instances we also used the numeri- 
cally exact quantum Monte Carlo method as an impu- 
rity solver, and generally found good qualitative and 
even quantitative agreement, supporting the validity of 
our IPT predictions in the relevant parameter ranges. 
Throughout the paper we express all energies in units of 
the bandwidth. 



III. PHASE DIAGRAM 

We first examine the evolution of the coexistence re- 
gion as disorder is introduced. Within this region, both 
metallic and insulating solutions are found, depending on 
the initial guess used in the iterative scheme for solving 
the self-consistency condition. Typical results are pre- 
sented in Fig. ^ showing the phase diagram obtained 
within DMFT-IPT at non-zero temperature, for varying 
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FIG. 1: Phase diagram for the disordered Hubbard model at 
non-zero temperature, (a) (U,T) diagram for different disor- 
der strengths, (b) (U,W) diagram at different temperatures. 
Ud and U C 2 lines are indicated in one of the plots, but similar 
definitions apply to the other results as well. 



levels of disorder W. For each level of disorder [shown in 
panel (a)] or temperature [shown in panel (b)], the first 
(from left) of the two lines, the so-called U&, indicates 
the stability boundary (i.e. the spinodal) of the insu- 
lating solution. Conversely, the second of the two lines, 
identified as U C 2, represents the boundary of the metallic 
solution. The coexistence region is found between these 
two lines, i.e. for U c \ < U < U C 2- Our results are in good 
quantitative agreement with previous results obtained in 
the T = limit in presence of disorder^ and also with 
non-zero temperature results in absence of disorder i As 
the disorder increases, the metal-insulator transition gen- 
erally moves to larger U. Physically, this reflects the fact 
that disorder broadens the bands and smears the gap, 
making it harder for the Mott-Hubbard gap to open, so 
that a larger U is necessary for the transition. At the 
same time, the temperature-dependent coexistence re- 
gion is found to shrink [Fig.^a)], persisting only below 
a critical end-point temperature T C (W) . At any given 
temperature, the principal effects of introducing disor- 
der [Fig. Bib)] are as follows: (1) both the U c i and U C 2 
lines move towards larger interaction potential; (2) the 
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lines become closer to each other as disorder increases. 
In fact, they both approach the W = U line as W — > oo. 

Having obtained these results in quantitative detail, 
we would like to understand the physical origin of this 
behavior. In the following we present simple analytical 
arguments relating the non-zero temperature aspects of 
the coexistence region to the evolution of its ground state 
properties. Our strategy is motivated by the following 
observations: (a) the shape of the non-zero temperature 
coexistence region [Fig.^a)] remains very similar at dif- 
ferent values of disorder; (b) its size, both in terms of 
temperature and in terms of [/-range, shrinks as disor- 
der increases. This suggests that the physical mechanism 
for the destruction of the coexistence region as the tem- 
perature increases is similar to that of the clean limit, 
where it is governed by decoherence processes due to in- 
elastic electron-electron scattering. Therefore, we begin 
our analysis by concentrating on the clean limit, where 
we show how simple estimates for the critical end-point 
temperature T c can be obtained. 



IV. COEXISTENCE REGION IN THE CLEAN 
LIMIT 



The coexistence region at non-zero temperature is de- 
limited by the two spinodal lines U c i (T) and U c2 (T) ; the 
critical end-point temperature T c is reached when these 
two boundaries intersect. To estimate T c using the T = 
properties of the model, we need to understand the tem- 
perature dependence of each of these lines. 



A. Insulating spinodal 

The insulating spinodal U c \ (T) essentially corresponds 
to the closing of the gap separating the two Hubbard 
bands in the Mott insulator. Its temperature depen- 
dence should thus reflect that of the Hubbard bands. 
In contrast to the correlated metallic state close to the 
Mott transition, the insulating solution is not character- 
ized by a small energy scale in the coexistence region. 
Accordingly, it is not expected to have strong temper- 
ature dependence; its weak temperature dependence re- 
flects activated processes across the Mott-Hubbard gap. 
Such activations only lead to (exponentially) weak round- 
ing/broadening of the Hubbard bands, which should very 
slowly reduce U C \{T) as temperature increases. Such be- 
havior is indeed clearly seen in our results. This tem- 
perature dependence is, however, much weaker than that 
characterizing U c2 (T). For purposes of roughly estimat- 
ing T c , to leading order we can ignore this weak temper- 
ature dependence, so that 



B. Metallic spinodal 

In the vicinity of the Mott transition, the metallic solu- 
tion is characterized by a low energy scale corresponding 
to the coherence temperature T* of a low-temperature 
Fermi liquid;* Above T* the heavy quasiparticles are de- 
stroyed, and the metallic solution becomes unstable. To 
estimate U c2 (T) we need to determine how this coherence 
temperature varies as the transition is approached. From 
detailed studies of the cleani and disordered^ Hubbard 
models within DMFT, it is known that this coherence 
temperature can be estimated as 



T* w AT F Z 



(5) 



where Tp is the Fermi temperature, A is a constant of 
order one, and Z is the quasiparticle (QP) weight defined 
as 



Z = 



1 - -ImEH 



(6) 



The behavior of Z is well known in the clean limit where 
it decreases linearly as U increases toward the metallic 
spinodal, viz. 



Z = C[U c2 (0) - U}. 



(7) 



From numerical studies^ the proportionality constant 
C « 0.45. Therefore, the coherence temperature can 
be written as 



T*(U) = ACT F [U c2 (Q)-U}. 



(8) 



We can now estimate the temperature dependence of 
U c2 (T) as that value of the interaction needed to set 
T*(U) = T, i.e. 

T = ACT F [U c2 (0) -U c2 (T)}. 

In other words 

U c2 (T) & U c2 (0) - BT, (9) 

where B = 1/ACTp. From our numerical results [see 
Fig. nja)] we find B w 22, giving A « 0.2, in reasonable 
agreement with estimates^ from the literature. 

Using these expressions for U C \{T) and U c2 {T), we ar- 
rive at the estimate for the critical end-point temperature 



T, 



[U c2 (0)~U cl (0)}/B, 



(10) 



U cl (T) « U c i(T = 0). 



(4) 



which agrees within 10% with our numerical results (see 
Fig.©. 



V. CRITICAL BEHAVIOR IN PRESENCE OF 
DISORDER 

Encouraged by the success of our analytical descrip- 
tion of the coexistence regime in the clean limit, we now 
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FIG. 2: (a) Average spectral function and (b) relative devi- 
ation of the distribution of pi(u), Ap/ (p), as a function of 
frequency for different values of the interaction potential. Ap 
is the standard deviation of the distribution of pi(uj), which 

where N is the 



is given by VEiCftM " <P*H»7(W ~ 
number of local site energies considered, 
used were T = 0.05 and W = 1.0. 
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Other parameters 



turn our attention to the effects of disorder. As in the 
clean limit, we would like to relate the finite temperature 
properties to the critical behavior of the quasiparticles at 
T = 0. To do this, we therefore concentrate on describing 
the critical behavior in presence of disorder. 

The principal new feature introduced by disorder 
within the DMFT scheme is the spatial variation of the 
spectral function, Pi(u>). This is shown in Fig. |2 at all 
energy scales: on the left we have the average spectral 
function and on the right the relative deviation of its dis- 
tribution, in the metallic phase. For each value of the 
interaction potential, the distribution of Pi(to) presents a 
large dip at lu ps and becomes broader as the frequency 
increases. This comes from the fact that at small fre- 
quencies the system is in the Fermi liquid regime. At 
finite temperature we observe the reminiscence of the 
perfect disorder screening seen at T = close to the 
Mott transition^ For large frequencies, the quasiparticle 
regime is no more valid and the appropriate description 
is in terms of Hubbard bands, resulting in an increase of 
the fluctuation in Pi(to). 

In the disordered case, the self-energy function 
presents site-to-site fluctuations, which lead to the spatial 
variations of the spectral function discussed above. The 
QP weights Zi = Z(si) now depend on the local site 
energy £<. To properly describe the approach to the Mott 
transition, we therefore must follow the evolution of the 
entire function Z(ei) as the transition is approached^ 




FIG. 3: (a) Quasiparticle weight as a function of the on-site 
energy for different values of the interaction potential as the 
U C 2 line is approached, for disorder strength W = f .0. The 
symbols are the numerical data, while the lines correspond to 
the fitting indicated in the plot, that is to a function with even 
exponents in e. (b) Fitted results for Z divided by Zo (the 
quasiparticle weight for e — 0) as a function of e, showing that 
close to the MIT the curves for different U scale. These results 
were obtained at a low but finite temperature T — 0.005. 



A. Behavior of local QP weights 

Given the self-consistent solution of our ensemble of 
impurity models, we calculate the local QP weights as 



Zi = 



1 - — ImEi(w) 

OOJ 



(11) 



Typical results are shown in Fig. Ota), where we plot 
Zi = Zfa) at T = 0.005, for disorder strength W = 1, 
as the metallic spinodal is approached by increasing the 
interaction U toward U C 2 « 1.9. We first observe that 
for small U, away from the transition, the QP weights 
Zi have strong dependence, with the smallest Zi at 
Ej = 0. Physically, this reflects the tendency for corre- 
lation effects (suppression of Z) to be the strongest on 
sites which are locally close to half-filling (singly occu- 
pied). Nonzero site energies favor the local occupation 
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FIG. 4: Distribution of quasiparticle weight for different val- 
ues of the interaction potential. The main plot shows how the 
curves collapse when we look at Z/Zq. The inset shows the 
results for Z itself. Other parameters as in Fig. [3] 

departing from half-filling, thus reducing the correlation 
effect, and increasing Zi. 

As U increases, all the Z»'s decrease, as in the clean 
case. But how does this affect the distribution of QP 
weights Zi — Z{ei)l At first glance it seems that the 
Si dependence becomes weaker, but a closer look reveals 
this not to be the case. As we shall now demonstrate, all 
the Zi's decrease linearly near the transition, i.e. they 
assume the form 

Z(U,e i )=K(e i )[U c2 -U}, (12) 

where only the prefactor K(ei) depends on £j. If, to 
leading order, these prefactors remain independent of the 
distance to the spinodal, then the entire family of curves 
Z{U,£i) can all be collapsed on a single scaling function. 
To verify this hypothesis, we define reduced QP weights 

Z*{ei) = Z(U,Ei)/Z(U,0). (13) 

If our scaling ansatz is valid, then the Z*(ei) should ap- 
proach a non-zero limit as U — ► U C 2, i.e. they should 
all collapse onto a single scaling function. As shown in 
Fig.[3Jb), this behavior is observed only for U sufficiently 
close to U C 2 [note that the data for U — 0.8 (further 
from the transition) show deviations from leading scal- 
ing]. This is precisely what we expect, since such simple 
scaling behavior typically occurs only within a critical 
region close to the metallic spinodal. 

B. Distribution P(Zi) of local QP weights 

Equivalently, we can characterize the QP weights by 
their probability distribution function P(Zi). Typical 
results for P(Zi) are shown in the inset of Fig. 0] As the 
Zi decrease near the transition, the distribution function 
P{Zi) changes its form and narrows down. However, if 



FIG. 5: Distribution of Z/Z for (a) smaller (W = 0.5) and 
(b) larger (W = 1.5) disorder than the one in 

Fig.H Other 

parameters used was T = 0.005. 



our scaling hypothesis is valid, then the shape of this 
distribution should approach a "fixed-point" form very 
close to the transition. More precisely, we expect the 
distribution for reduced QP weights P(Z*) to collapse to 
a single scaling function close to t/ c2 . Results confirming 
precisely such behavior are presented in Fig. 01 

An interesting question relates to the precise form of 
the fixed-point distribution function P(Z*), and how it 
may depend on disorder. In the clean limit, obviously, 
it reduces to S(Z* — 1) indicating that spatial fluctua- 
tions are suppressed. As the disorder increases, P{Z*) 
becomes very broad (as shown in Fig. |SJ, reflecting large 
site-to-site fluctuations in the local QP weights. This 
behavior may be regarded as a precursor of electronic 
Griffiths phases^ which emerge for stronger disorder, as 
found within statDMFT approaches ^ 

In essential contrast to the clean limit, the approach 
to the Mott transition in presence of disorder thus needs 
to be characterized by the entire probability distribution 
function of QP parameters. At first glance, this may 
appear to require a description considerably more com- 
plex than in the absence of disorder. However, we have 
demonstrated that in the critical region the distributions 
approach a fixed point form, allowing for "single parame- 
ter scaling," in close analogy to the clean Mott transition. 
This finding immediately suggests that our arguments 
describing the finite temperature coexistence behavior in 
the clean limit may successfully be extended to the dis- 
ordered case as well, allowing for a complete qualitative 
description, which we discuss in the following section. 

VI. COEXISTENCE REGION IN PRESENCE OF 
DISORDER 

Within the DMFT formulation, the disorder is not 
expected to qualitatively affect the temperature depen- 
dence of the insulating spinodal, since the forms of the 
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FIG. 6: Average quasiparticle weight as a function of the 
interaction potential for different values of disorder. 



Hubbard bands remain qualitatively similar to that in 
the clean limit. The principal effect of disorder in the 
Mott insulating phase is to simply broaden the Hubbard 
bands, which retain well defined (sharp) band edges due 
to the CPA-like treatment of randomness in the DMFT 
limit. Indeed, our quantitative results [see Fig.^a)] con- 
firm that Ud (T) ss U c i (0) retains very weak temperature 
dependence, as in the clean case. The only modification 
is that E/ c i(0) rapidly grows as disorder is increase, re- 
flecting the disorder-induced broadening of the Hubbard 
bands. 

The metallic solution is again found to be unstable 
above a certain coherence temperature T*(W, U), which 
defines the locus of the metallic spinodal U C 2(T). An 
added subtlety is that different sites start to decohere at 
different temperatures, an effect that earlier work 9 found 
responsible for a nearly-linear temperature dependence 
of the resistivity in the disordered metallic phase. Nev- 
ertheless, sufficiently close to the Mott transition (within 
the coexistence region), a sharply defined temperature 
scale T* (W, U) emerges where the metallic solution sud- 
denly disappears and where the qualitative form of the 
spectrum changes on all sites. This temperature scale 
defines the locus of the metallic spinodal, corresponding 
to the equation 

T = T*(W,U c2 ). (14) 

At first glance, it is anything but obvious how 
T*(W, U) should be estimated. As in the clean case, the 
reduction of this temperature scale as the transition is 
approached must reflect the behavior of the local quasi- 
particle weights Zi, and presumably depend on the pre- 
cise form of the distribution function P(Zi). As we have 
seen, however, all the local QP weight scale in a simi- 
lar fashion in the critical regime, which suggests that a 
reasonable estimate may be obtained simply from their 



FIG. 7: Temperature at which the U c i and U C 2 lines merge in 
the (U, T) phase diagram as a function of disorder. The plot 
shows both the results obtained directly from the numerical 
data as well as those calculated from the linear fitting to the 
U C 2(T) line. In the latter, T c was calculated using the values 
of Ud at T — 0.005, except for W = where the result at 
T = 0.0075 was used. 

average value 

(Zt) = J de l P{e l )Z l . (15) 

At least for sufficiently weak disorder, we may expect 
that [cf. Eq. (5)] 

T*(W, U) ~ ATp {Zi) , (16) 

where A « 0.2 as in the clean case. Using the fact that 
all ZiS decrease linearly near the transition, we expect 

(Zi) = C(W)[U c2 (0) - U]. (17) 

To confirm this, we explicitly calculated (Zi) as a func- 
tion of U for different levels of disorder; the results are 
shown in Fig. El We conclude that C(W) w C(0) « 0.45. 
These results suggest that the metallic spinodal should 
take the form, 

U c2 (W, T) « U c2 (W, 0) - B(W)T, (18) 

where B(W) « B(0) = 22. Our non-zero temperature 
results for U c2 (W, T) [see Fig. EI a)] fully confirm these 
expectations. Based on these results, we finally obtain 
the desired expression for T C (W) of the form 

T C (W) « [U c2 (W, 0) - U cl (W, 0)]/B(0). (19) 

To test the proposed procedure, we have used the values 
for U C \(W) and U c2 (W) at the lowest temperature of our 
calculation (T = 0.005) to estimate T C (W). As we can 
see from Fig. our analytical estimates are found to 
be in excellent agreement with results of explicit non- 
zero temperature calculations. The decrease of T C (W) 
with disorder thus directly reflects the "shrinking" of the 
coexistence region at low temperature, which in its turn 
reflects the decrease of the energy difference between the 
metallic and the insulating solution. 
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VII. CONCLUSIONS 

In this paper we have used a DMFT approach to ex- 
amine the effects of disorder on the critical behavior near 
the Mott metal-insulator transition, with special empha- 
sis on non-zero temperature properties associated with 
the two spinodal lines U c \ and U C 2- By using a combi- 
nation of numerical results and analytical arguments we 
have demonstrated that simple scaling behavior emerges, 
providing a complete description of the critical regime. 

In contrast to the clean case, the presence of disorder 
requires one to examine the entire distribution of local 
spectral functions, Pi(w), describing how the local spectra 
varies with position in the sample. This can be probed 
with scanning tunneling microscopy (STM). Notice that 
the distribution function describing the site dependence 
of Pi{uj) will depend on the frequency of observation: it 
will be broader at higher energies [as seen in Fig. Hfb)], 
where a real space picture is appropriate to describe the 
Hubbard bands, and narrower at low frequencies, where 
a quasiparticle description in k space is appropriate. This 
is a manifestation of frequency dependence of the disor- 
der screening discussed in an earlier paper by some of 
us£ 

In the metallic regime, at low temperatures, the spec- 
tral function can be parametrized in terms of the distri- 
bution of quasiparticle parameters, which displays sim- 



ple scaling properties. This allowed us to characterize 
the behavior near U C 2 using a single parameter scaling 
procedure. The approach to U C 2 thus retains a character 
qualitatively independent of the level of disorder, where 
the vanishing of quasiparticle weight signals the trans- 
mutation of itinerant electrons into localized magnetic 
moments. 

Within the examined DMFT formulation, the region 
between the two spinodal lines U& and U C 2 although re- 
duced in size and extent cannot be completely eliminated 
no matter how large the disorder. Of course, these predic- 
tions are applicable only for weak enough disorder where 
Anderson localization effects can be ignored. Extensions 
of DMFT that incorporate Anderson localization mecha- 
nisms at zero temperature are available jifiiil but applying 
these approaches to examine the non-zero temperature 
behavior near Mott- Anderson transitions remains an in- 
teresting research direction. The behavior at the first 
order transition line and the actual nucleation of cither 
the metallic or insulating phase, between U c i and U C 2, 
are also strongly modified by disorder, and this as well is 
left for future study. 

The authors thank A. Georges and D. Tanaskovic 
for useful discussions. This work was supported by 
NSF grants DMR-9974311 and DMR-0234215 (VD) and 
DMR-0096462 (GK). 
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